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Abstract. - We study numerically and partly diagrammatically the renormalization of the 
electron-spin interaction or vertex in a two-dimensional one-band Hubbard model with spin- 
fluctuation momentum transfer q = (7r,7r). We find that the renormalized electron-spin vertex 
decreases quite generally with decreasing temperature at all doping densities. As a combination 
of two concurring effects, i.e. the decrease of the vertex and the increase of the spin suscepti- 
bility, the effective pairing interaction increases with lowering temperature in the intermediate- 
correlation regime, but flattens off in the strong-correlation regime. Our findings indicate that 
in the high-T c cuprates the pairing mediated by antiferromagnetic spin fluctuations is substan- 
tially reduced due to vertex corrections. 



Introduction. - The spin-fluctuation mediated interaction between charge carriers is a 
corner stone in a variety of theoretical and experimental issues [1,2], which are presently in 
the center of focus of the high-T c superconductivity research: it has been proposed that a 
spin-1 resonance mode, which is prominent in the magnetic response measured by neutron 
scattering, is responsible for significant features ( "kink" ) seen in other spectroscopies such as 
photoemission [3-5] and optical absorption, which are sensitive to the charge dynamics, and 
even that this mode is the boson responsible for "mediating" the superconducting pairing [6] . 
This also has led to a variety of counter-arguments based on the small spectral weight of the 
resonance mode and (assumed weak) coupling to electron-hole pairs [7]. However, neither the 
original proposals for the resonant spin-fluctuation mode scenarios nor the counter-arguments 
take into account the effects of vertex corrections. The vertex function 7(k, q) = 3kq/.9kq' 
which gives the ratio of the coupling <?kq of the boson (q = (tt, it) magnon or resonance 
mode) with a dressed quasiparticlc (which includes vertex corrections and also quasiparticle 
renormalization) to the bare vertex , is expected to be substantially renormalized by 
strong electronic correlations. This has recently been demonstrated for the electron-phonon 
(el-ph) vertex [8] by using numerically accurate quantum-Monte-Carlo (QMC) techniques [9]: 
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whereas at weak and intermediate Coulomb (Hubbard U) interactions, screening is the dom- 
inant correlation effect suppressing the el-ph coupling, at larger U values irreducible vertex 
corrections become dominant and give rise to an unexpected and substantial increase in the 
forward direction. 

The need for clarifying the role of vertex corrections in the high-T c schemes involving the 
exchange of spin-fluctuations also derives from a rather dramatic finding by Schricffcr [10]: 
he argued that the spin-fluctuation vertex is strongly suppressed in the long-ranged anti- 
ferromagnetic (AF) state or, in the paramagnetic state but with strong AF precursors, i.e. 

7 (k, q) <x [(q - Q) 2 + jzl* ~ . - 1 -7T , (1) 
e [ X (q,c; = 0)]2 

where £ stands for the AF correlation length and %(q, u> = 0) is the spin susceptibility. This 
finding is further confirmed by Chubukov et al. [11]. The vanishing of 7 at q = Q = (7r,7r) in 
the long-ranged AF state is a consequence of the Adler principle: in the (AF) ordered state, 
the magnons are Goldstone bosons and, therefore, they always decouple from other degrees 
of freedom [11]. This then leads to a strong renormalization of the pairing interaction V in 
the single spin-fluctuation exchange approximation: 

y(k,q)«| 7 (k,q)| 2 X (q), (2) 

which instead of diverging tends to a finite positive constant as q approaches Q and £ — > 00. On 
the other hand, for the appearance of d-wave superconductivity one would need a substantial 
enhancement of the interaction near Q. There remains the question of how large the AF 
correlation length £ must be for this reduction effect to be substantial. Furthermore, the 
theoretical argument uses the "frozen in" electron-spin model, where the holes move in an 
unaffected (by the presence of holes) AF background. 

Using QMC techniques [9] , we aim at clarifying the role of electronic correlations on the 
effective electron-spin fluctuation (el-sp) coupling in a two-dimensional one-band Hubbard 
model. By comparing QMC results at a typical low doping (~ 10%) with diagrammatic 
(~ U 2 ) calculations, we infer that the "weak correlation" regime, i.e. where 7 ~ 1, is confined 
to rather small U (U < 2t, t: hopping) values (see Fig. [2J. In this regime, and also in 
the intermediate correlation regime (U ~ it) where the bare vertex is already substantially 
suppressed (~ 50%), the strong increase of the susceptibility x with decreasing temperature 
dominates the effective pairing potential V in Eq. Q. Thus, V(p,q = (7r,7r)) with p close 
to the Fermi surface increases for decreasing temperatures. However, this standard behavior 
even qualitatively changes in the physically relevant strong correlation (U > St) regime: 
here, we do find in the underdoped regime that the vertex 7 is so substantially reduced at low 
enough temperatures (T < J = At 1 /U) that this reduction over-compensates the simultaneous 
increase of the spin susceptibility. Thus, at temperatures below the magnetic scale J, the 
vertex reduction introduces a substantial suppression of the effective pairing potential. At 
a first glance, this finding displays a similarity to Schricffer's observation (see Eq. |5J). We 
will argue, however, that in our calculation the AF precursor has not developed and that 
our findings are not consistent with a picture of hole fermions embedded in a long-range AF 
background, but rather with a short-range "real-space" picture, valid for strong correlations. 
Here, the spin fluctuation couples to the quasiparticlc which forms a "spin-bag" , i.e. a hole 
dressed by short-range AF fluctuations. 

Model and numerical approach. - Our starting point is the one-band Hubbard model, 

(ij) ,<r i 
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The operators c| CT and c itT as usual create and destroy an electron with spin er at site i, 
respectively and the sum (ij) is over nearest-neighbor lattice sites. Here, U is the onsite 
Coulomb interaction and we will choose the nearest-neighbor hopping t as the unit of energy. 

In our simulations, we have used the linear-response technique as in Ref. [8] in order 
to extract the el-sp vertex. In this method, one formally adds to Eq. © the interaction 
with a momentum- and (imaginary) time-dependent spin-fluctuation field in the z-direction 
Sqe- iqoT [12] in the form [13] 

Hel-s P = £ <?V4 +q(r C k<7 S « e " ?9 ° T > ( 4 ) 
kqcr 

where g£ q is the bare el-sp coupling (equal to the Hubbard U in the one-band Hubbard 
model). One then considers the "anomalous" single-particle propagator in the presence of 
this perturbation defined as [12] 

G A (p,q) = -fdr e^ + ^^T T c p+ ^(T)cl a (0))H + H el _ ap , (5) 
Jo 

Here () H+H ci - sp is the Green's function evaluated with the Hamiltonian H + H e i_ sp . Dia- 
grammatically, Ga(p, q) has the structure shown in Fig. ^so that the el-sp vertex T(p, q) can 
be expressed quite generally in terms of Ga and of the single-particle Green's function G(p) 
in the form 

r( \-v 1 1 1 G A ( P ,g) 

1 S%gl q S (l l + U Xzz ( q )G{p + q)G(p) ' W 

with Xzz (q) the longitudinal spin susceptibility. Due to the limit procedure in Eq. it is 
sufficient to calculate the leading linear response of Ga to i? e i-sp, which is given by 

G A ( P ,q) = 5 q f dre***** f dre^' £ ff ° q x 
Jo Jo kCT , 

(Ttv'cI +c1 At' + + ) Ck<T ,(^)c P+q .(r)ct >CT (0)) if , (7) 

where + is a positive infinitesimal. Notice that 5 q cancels in Eq. El The two-particle Green's 
function in Eq. (J7J) is evaluated with respect to the pure Hubbard Hamiltonian (Eq. @ ) . 

The effective vertex 7 describing scattering processes between quasiparticles and spin waves 
is obtained by taking into account the wave-function renormalisation Z(p): 

j(p,q)= , TiP 'f . (8) 

Numerically, Z is obtained as Z(p) = lm[l/G(p)]/po [8,12]. In terms of 7(^,9) and of the 
spin susceptibility \ zz (q), the effective pairing interaction V under the exchange of a single 
spin wave is expressed in the form 

V(p,q) - h( P ,q)\ 2 U 2 Xzz (q) 

= [|Re7(p,9)| 2 + |Ini7(p,g)| 2 ]C/ 2 X ^('?), (9) 

with 

1 



= \ f dr 



and 
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G A (p,q)= 




Fig. 1 - Diagrammatic representation of GA(p,q) within linear response to S q . The thick solid lines 
represent dressed single-particle Green's functions of the Hubbard model. The wavy line denotes the 
external perturbation in Eq. (J1J. The dashed line represent the Hubbard interaction U and the closed 
ellipse stands for the longitudinal spin susceptibility Xzz(g). 



Cko 



(10) 



where V(p, q) contains the contributions from both the real and imaginary parts of the vertex 
7. Notice that due to the renormalisation in Eq. [SI V(p, q) describes the effective interaction 
between quasiparticles. 

Our numerical Monte Carlo simulations were performed on an 8 x 8 lattice at different 
doping densities and different temperatures [12]. In the high-T c cuprates, the charge carriers 
near the (n, 0) region are strongly affected by antiferromagnetic spin fluctuations. Therefore, 
we will examine the particular scattering process in which the incoming electron and spin 
fluctuation carry momenta p = (—it, 0) and q = (n, 7r), respectively. Within our p-points 
mesh, these points lie sufficiently close to the Fermi surface. 

Results. - The diagrams contributing to the vertex T up to order U 2 are displayed in 
Fig. Ela). Figs. EJb) and^c) show the comparison of QMC calculations with perturbation 
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Fig. 2 - (a) Low-order Feynman diagrams for the el-sp vertex T(jp,q). The thin solid lines are the 
non-interacting Green's functions and the dashed lines represent the Hubbard interaction U. The 
wavy lines stand for the spin-fluctuation fields. Real (b) and imaginary (c) parts of F(p, g) vs U 
at 5 = 0.12 and f3 — 2. The solid circles are Monte Carlo results and the open symbols show the 
perturbation theory contributions shown in [51 a). 
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U=4; 5=0.12 (circles) 
U=8; 5=0.12 (squares) 
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Fig. 3 - Real (a) and imaginary (b) parts of the effective el-sp vertex 7(p, q) vs T at U = 4 and U = 8 
for the doping density S = 0.12. Figs, (c) and (d) show the T-dependence of the spin susceptibility 
Xzz{q) and the effective pairing interactions between quasiparticles V(p, q). 



theory. As one can see, the perturbative calculations are in good agreement with Monte Carlo 
simulations only for U < 2. 

Monte Carlo results for j(p,q), Xzz(q), and V(p,q) are displayed in Figs. Et a )"Etd). We 
notice that in the undcrdoped regime (S = 0.12) for both intermediate correlation (U = 4) and 
strong correlation (U = 8), 7 is strongly rcnormalized at characteristic temperatures (below 
T « J = 0.5 for U = 8). There is a clear temperature trend observable in Re7, which tends 
to go to zero (or to a very small value), at least for U = 8. At our lowest accessible T's, 
Im7(p, q) is small at weak and intermediate correlations, but for U = 8 it can become larger. 
Our numerical results clearly show that, whereas in the weak and intermediate correlation 
regimes the strong enhancement of the spin susceptibility dominates the low-T behavior of 
V, which then strongly increases for decreasing T, the behavior is different in the strong 
correlation (U = 8) regime. Here, the pairing potential V only displays a mild increase or 
a saturation at low T's. Nevertheless, even at the lowest temperatures accessible at U = 8, 
our conclusion that vertex corrections are crucial is already quite clear: there is an order of 
magnitude reduction [~ O(10)] in the effective pairing interaction V (Fig. Old)) compared 
with the RPA (with full \zz) result. Our finding at U = 4 is in good agreement with the work 
of Bulut et al. [14] , which shows that the effective particle-particle interaction created by the 
Hubbard U increases with lowering temperature and can reach large values. 

Figure 0] gives for both the intermediate (U = 4) and strong correlation (U = 8) cases the 
doping dependence of the susceptibility Xzz(l), of the renormalization factor \j(p, q)\ 2 , and 
of the effective pairing interaction V(p,q). Again, we note the competing influences of the 
temperature on \ an d 7: X increases both as a function of lowering doping and temperature; 
this happens in both U = 4 and U — 8 cases. On the other hand, the effective el-sp vertex 
7 decreases as a function of lowering temperature (as we have seen in Fig. EJ) , and is almost 
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Fig. 4 - Spin susceptibility Xzz(?)i the renormalization factor |7(p,g)| 2 (i. e. the reduction of V 
with respect to the RPA result), and the effective pairing interaction V(p, q) as a function of doping 
density S for (a) U = 4 and (b) {7 = 8. The value of the inverse temperature /3 is indicated by the 
shape of the symbol. 



doping independent. As already mentioned, at intermediate correlations, the T-dependence 
of x still dominates the pairing interaction, resulting in an increase of the effective pairing 
upon reducing T, especially for low dopings. But even here, at U = 4, vertex corrections 
reduce V substantially, as can be seen in the \-f\ 2 value, which is reduced compared from the 
bare value 1 down to about ~ 0.25. This vertex influence on the effective pairing interaction 
is most dramatic at U = 8. Here \j\ 2 reduces V by a factor of 8 at the lowest accessible T's 
and dopings 5. 

Discussion and conclusions. - How can we understand the rather dramatic role of vertex 
corrections in the strongly-correlated underdoped regime, where we have a short-ranged (of 
order Cu-Cu distance) correlation length £ [15] and not Schrieffer's situation, where holes 
move in an unaffected (by the charge carriers) AF (£ — > oo) background? The important 
point is that in our calculation, a strong vertex suppression is obtained even though the AF 
precursor has not developed. What our QMC calculations demonstrate, is that below the 
characteristic temperature T ~ J, vertex corrections become so large that in the physically 
relevant strong-correlation (U = 8) regime, an order of magnitude reduction in the pairing 
interaction results. This reduction happens more or less independent of doping. This doping 
independence is at variance with Schrieffer's result where of course £ depends crucially 
on doping. Thus, our results suggest the following picture: It is known from our QMC 
work on the single-particle spectral function A(k,uj) [15] that below T ~ J, a "band" of 
width O(J) forms, where "spin-bag" -like quasiparticles propagate coherently. This happens 
again more or less independent of doping, i.e. from the underdoped insulator (£ — > oo) 
up to optimal doping (£ wCu-Cu distance). It can also be shown on the basis of exact 
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diagonalization that it is the "spin-bag" quasiparticle and not the incoherent "lower Hubbard 
band" background which couples most effectively to the perturbing spin potential H e i- sp in 
Eq. (0J [16]. Therefore, as soon as the quasiparticle with its spin dressing has been formed 
(for low enough T"s and large enough J7's) the scattering amplitude or vertex is more or 
less independent of doping (see Fig. This picture is clearly different from the doping 
dependent "screening effect" implicit in Eq. Q . Our results were obtained for a model with 
nearest-neighbor hopping only. For a nonvanishing next- nearest-neighbor hopping t' jt < 0, 
spin fluctuations get damped. In this case, calculations within the spin-fermion model show 
that vertex corrections are reduced [11]. We carried our QMC calculations for finite t' jt = 
—0.25, showing that indeed vertex corrections arc smaller than for t'/t — [17], although the 
difference is not as pronounced as in Ref. [11]. 

In summary, based on quantum Monte Carlo simulations, we have studied the renormal- 
ization of the cl-sp interaction or vertex in the one-band Hubbard model. We found that 
the renormalized el-sp vertex decreases quite generally with decreasing temperature. On the 
other hand, the temperature dependence of the effective pairing interaction is rather different 
in the intermediate- and strong-correlation regimes: It increases with lowering temperature 
in the intermediate-correlation regime, but flattens off in the strong-correlation regime. 
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